Electronic compressibility of graphene: The case of vanishing electron correlations 

and the role of chirality 
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A recent surprising finding that electronic compressibility measured experimentally in monolayer 
graphene can be described solely in terms of the kinetic energy [J. Martin, et al., Nat. Phys. 
4, 144 (2008)] is explained theoretically as a direct consequence of the linear energy dispersion 
and the chirality of massless Dirac electrons. For bilayer graphene we show that contributions to 
the compressibility from the electron correlations are restored. We attribute the difference to the 
respective momentum dependence of the low-energy band structures of the two materials. 
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In an interacting electron system of uniform density, 
the (inverse) electronic compressibility oc dfi/dn 
(where \i is the chemical potential and n is the elec- 
tron density) is a fundamental physical quantity that 
is intimately related to the strength of interelectron 
interactions 1,2 . First measured in 1992, the compressibil- 
ity of a two-dimensional electron gas 3 provides valuable 
information about the nature of the interacting ground 
state, particularly in the strong-coupling regime where 
(in addition to the exchange energy) the Coulomb in- 
teraction energy is also known to play a dominant role. 
In this context, a recent report on the measurement 
of electronic compressibility in monolayer graphene 4 re- 
vealed behavior which was totally unexpected 5 . In this 
work, scanning single-electron transistor microscopy was 
used to measure the change of local electrostatic po- 
tential (and thereby change in local chemical poten- 
tial) of a graphene sample when the carrier density was 
modulated 5 . The observed results for the local inverse 
compressibility were found to be quantitatively described 
by the kinetic energy alone and the authors speculated 
that the exchange and correlation energy contributions 
to the compressibility either each other cancel out or are 
negligibly small. This interesting puzzle has remained 
unsolved because the approximate theoretical schemes 
adopted by various authors to investigate electron corre- 
lations in graphene 6 ' do not find any such cancellations. 
Similarly, the recently reported Hartree-Fock studies of 
compressibility 8 in monolayer and bilayer graphene do 
not consider electron correlations at all and are therefore 
not in a position to address this important issue. 

In this paper, we investigate the role of electron corre- 
lations in monolayer and bilayer graphene. We show how 
in monolayer graphene, two fundamental properties of 
the system, viz., the linear energy dispersion and chirality 
conspire to allow the exchange and correlation contribu- 
tions to vanish, just as was observed in the experiment 5 . 
In bilayer graphene on the other hand, where the low- 
energy quasiparticles are massive chiral fermions 9 ' 10 , the 
parabolic dispersion does not allow this vanishing of the 



two energies, and the kinetic energy retains a dependence 
on the electron correlation function which manifests in 
the electron compressibility. 

The low-energy charge carriers in monolayer graphene 
behave as massless Dirac fermions described by the 
single-particle Hamiltonian 7^ m oc <r • p 1; which is linear 
in momentum p and where the subscript refers to the co- 
ordinate label of the electron on which it acts. The eigen- 
states of the Hamiltonian are uniquely labeled by quan- 
tum numbers representing the wave vector q = p/h, the 
band (conduction/ valence) b, valley (pseudospin) £, and 
the z-component of the real electron spin ex. The wave 
functions are of the form, ^(r) = e ,q r x, where x is an 
eight component spinor 11 . A full analytical study of the 
many-electron system in graphene is clearly an impossi- 
ble task. However, most of the clues to the puzzle involv- 
ing the measured compressibility described above can be 
found at the level of two electrons, which is amenable to 
a fully analytic solution. We therefore start with a two- 
electron system where the electrons occupy the states a 
and (3 corresponding to the full sets of quantum numbers 
(cia,/3> b at /3, ^ a ,/3j °a,/9)- We denote by ip the antisymmet- 
ric noninteracting two-electron wave function 

P(ri,r 2 ) = -j= [il> a (ri)ipp(r 2 ) - ^g(ri)V>a(r 2 )] • 

The correlations due to the mutual Coulomb interaction 
are introduced by multiplying the free-particle wave func- 
tion by a generic correlation factor F as 

# = F(ri,r 2 V(ri,r 2 ). 

At this stage, a precise definition of F is not neces- 
sary. The only requirements are that it should be a real 
function, and to preserve the antisymmtery of the cor- 
related wave function '5 it is assumed to be symmetric 
with respect to exchange of the particle coordinates, i.e., 
F(r 1 ,r 2 ) = F(r 2) r 1 ). 
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In order to evaluate the two-particle energy we have to normalize the wave function 'J. A straightforward calculation 
yields 

ll*H 2 H*l*H / <fridr 2 F(ri,r 2 ) 2 {l-^ o£/J [l + b a b fj cos{9 a - 9 fj )\ cosQ}, 

where Q = (q^ — q a ) • (ri — r^), and a ^ are the polar angles of the momenta q a ,/3- Evaluation of the expectation 
value of the kinetic part of the Hamiltonian T m = {^[T™ + T 2 ra \^) leads to the expression 



= \hv F J dv 1 dviF^iS^S^ [eVzip^l + babpe-W'-oj) + Z} a (l + b a b e^ e ^) 

+ e-^Z\ a (l + babpeWe-'J) + e~^Z% (l + bjo^-^) 
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where F stands for F(ri,r 2 ), v F is the Fermi velocity, and Z 3 af} is shorthand for 



dF 



Due to the linearity of T m in the momentum operators, 
only the first-order derivatives appear in the integrand. 
Terms in T m of the form 



dr—F 2 

3 dx/ ' 



clearly vanish due to the antisymmetry of the integrand. 
Most of the terms left after the volume integration can- 
cel each other as a consequence of the spinor structure 
of the single-particle wave functions, which is a direct 
manifestation of the chirality of the electrons. The only 
surviving terms sum to 

T™ = hv F (b a q a + b [m )m\\ 

that is, the kinetic energy expectation value (T m ) is sim- 
ply the sum of the single free-particle kinetic energies, 
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Tiv F (b a q a + b p q p ) = T m 
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and does not depend on the correlation function F at all. 
We expect similar cancellations for higher electron num- 
bers, although analytical expressions become intractable 
already at the level of three electrons. 

Complete cancellation of correlation contributions to 
the kinetic energy (never observed in conventional elec- 
tron systems) creates an unusual situation as we shall 
now describe. In the thermodynamic limit, the potential 
energy (per particle) V is usually expressed in the form 

{V)=nJ dv[g(r)-l}V c<m :(r), 



where n is the single particle number density, Vboui is 
the Coulomb potential and g(r) is the pair-correlation 
function which, for r = ri — r 2l is given by 

9(|r|) = w x 



x J dr 3 . ..dr N |*(ri,r 2 ,r 3 , . . . , rjv)j S 



where N is the total number of electrons. The energy 
functional (per particle) £ m is now 



dr[g(r)-l]V Coul (r), 



where = T™/N is the kinetic energy per particle. Its 
variation with respect to g(r), an essential step in deter- 
mining the optimal g(r), would yield an unusual Eulcr- 
Lagrange equation, VcouiM = 0, which is clearly not 
the case in graphene 12 . To resolve this dilemma we note 
that the energy functional £ m is actually not bounded 
below: we can choose correlations such that the poten- 
tial energy takes arbitrarily large negative values. This 
implies that to determine the optimal g(r) the energy 
functional derived above is not sufficient and additional 
physical constraints, for example, that g(r) should cor- 
respond to the correct number of states in each band, 
would be necessary. Clearly, determination of the opti- 
mal pair-correlation functions for massless Dirac fermions 
in graphene is a nontrivial problem 13 . However, we be- 
lieve that the expression for the functional £ m is of the 
correct form; i.e., once the correct pair-distribution func- 
tion g{r) is found, one could evaluate the correct energy 
from the above form of the energy functional. 
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Let us now turn our attention to the compressibility as 
defined in the introduction. To that end we first evaluate 
the variation 8 n £ m of £ m with respect to n: 

S n £ m = ^Sn + SnJdv [g(r) - l]V Co ui(r) 

+ n [ drVcoui^r^-Sn. 



Sn 



From this we can read the derivative as 

d£ m dt f J . . . ilT „ j' 5g{r) 

on on J J On 

The compressibility will then be proportional to 
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+ n / rfrVboui 



S(6g(r)/Sn) 
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We therefore need to make an assumption or estima- 
tion of the functional derivative Sg(r)/5n. The pair- 
correlation function is a many-body quantity, so its exact 
evaluation is impossible. Also, it is not possible to calcu- 
late functional derivatives numerically, and any analyti- 
cal approximation will necessarily obscure the true func- 
tional dependence that we require. In conventional two- 
dimensional electron systems, g(r) varies only slightly as 
a function of density except at very low densities where 
it starts to develop a prominent peak as a precursor to 
Wigner crystallization 14 . Therefore, we appeal to pre- 
vious work which shows that there are no phase tran- 
sitions (such as Wigner crystallization) as the density of 



the graphene system is varied 15 and so we expect that the 
functional variation of g(r) with the density will be neg- 
ligable in this system. Alternatively, we could consider 
a slightly less stringent condition J dr Vboui Sg ^ = 
which implies that the interaction energy depends lin- 
early on the density of Dirac electrons. Under cither of 
these assumptions, the compressibility is described en- 
tirely by the kinetic energy 



d 2 £ m 
dn 2 



dn 2 ' 



in accordance with the experimental observation 5 . In 
arriving at this striking result, there are two basic prop- 
erties of monolayer graphene that play crucial roles: the 
linear energy dispersion and chirality of massless Dirac 
electrons. 

This immediately invites the question: What happens 
in bilaycr graphene 9,10 , where the low-energy charge car- 
riers behave as massive chiral fermions and as such the 
Hamiltonian is quadratic in momentum operators near 
the charge neutrality point. The single-particle Hamil- 
tonian is Ti = — 2m T ( cr ' Pi) <T a:(cr ■ Pi) where m* is the 
effective electron mass generated by the inter layer cou- 
pling, and has the spectrum e a = (% Q ) 2 associated with 
it. Evaluating the expectation value of an arbitrary two- 
particle wave function as in the monolayer case, an inter- 
mediate expression for the kinetic part of the two-particle 
energy is 



dF dF\ 
cos(20 Q - 6 f3 )] — + [sin Bp + sin(20 Q - 6 f3 )} —j 



T h = i(e Q + ep) - J dridr 2 ^^[iqp ([cos 6 

(dF dF \ ~i 

[cos 6 a + cos(9 a - 29,3)} + [sm9 a - sm(9 a - 26p)] — j j 

/h 2 F r ( dF dF\ 

dr x dr 2 —e lQ {iq a ( [cos(9 Q + cos{6 a - 26p)] — + [sin6» Q - sin(6» Q - 29 p )} —J 

(dF OF \ 1 

[cos Bp + cos(26» a - Bp)] — + [sinOp + sin(2(9 Q - Op)] -^-) } 



where we have already excluded terms containing mixed 
second derivatives of F which are identically zero on in- 
tegration, and those which trivially sum to zero. The 
integrals of terms with single derivatives of F are finite, 
and the prefactors (coming from the pscudospinor part of 
the products of wave functions in the expectation value) 
do not cancel each other as they did in the monolayer 
case, but constructively sum to a finite result. This non- 



cancellation is a feature of the sublattice structure of the 
electronic wave function in bilayer graphene, resulting 
from the quadratic nature of the low-energy dispersion 
relation. On evaluation of the remaining integrals, and 
after some elementary algebra, the energy functional is 
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found to be 

X [ill + <$) cos (^ - dp) - 2q a qp] 

where F is the Fourier transform of the correlation func- 
tion. We can easily see that there is a nonzero contribu- 
tion from the electron correlations to the kinetic energy 
in this functional, and therefore taking the derivatives 
with respect to n yields a compressibility which depends 
nontrivially on them. It is clear that this additional term 
will also be present in the many body energy, as it's inte- 
gral over momentum is manifestly finite. We would also 
expect that for bilayer graphenc where the excess elec- 
tron density is high enough that the Fermi energy is in 
the energy range where the linearity of the spectrum is 
restored, that the effect of the correlations in the energy 
functional will again be supressed. Quantitative com- 
putation of this term requires precise knowledge of the 
radial dependence of F, and the relation between F(r) 
and g(r) 16 . Both of these issues are beyond the scope of 
the present paper, however, as indicated in our present 
results for bilayer graphene, experimental observation of 
a shift in compressibility from the pure kinetic energy 



contribution would provide a way to directly determine 
the strength of electron correlations in that system. 

To conclude, we have demonstrated that in monolayer 
graphenc, the electron correlations analytically vanish 
from the two particle kinetic energy. This, and the neg- 
ligability of Sg(r)/Sn lead to the absence of the elec- 
tron correlation function in the compressibility, as seen 
in recent experiments 5 . Conversely, the restoration of a 
quadratic band structure in bilayer graphene means that 
the correlations are present in the kinetic energy func- 
tional and compressibility in this case. Our work strongly 
suggests that quantitative agreement between the single 
particle theory and the experimental results in mono- 
layer graphene has its origin in the fundamental proper- 
ties (and in particular, the linear band structure) of this 
system. Also, the experimental system exhibits nonho- 
mogeneity of the charge distribution, so a full many-body 
calculation of the compressibility would have to include 
this detail. However, the cancellation of the correlation 
function from the monolayer two particle energy func- 
tional is independent of the energy of the electrons, and 
will therefore persist in the inhomogeneous system. 
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